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SYMBOLIC COMPUTATIONS IN DIFFERENTIAL GEOMETRY 



DIEGO CONTI 



00 

^~^ Abstract. We introduce the C++ library Wedge, based on GiNaC, for sym- 

^^ bolic computations in differential geometry. We show how Wedge makes it 

^vj possible to use the language C++ to perform such computations, and illus- 

trate some advantages of this approach with explicit examples. In particular, 
we describe a short program to determine whether a given linear exterior dif- 
ferential system is involutive. 



o 
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Introduction 



o 

\—A There are many computationally intensive problems in differential and Riemann- 

r^ ian geometry that are best solved by the use of a computer. Due to the nature of 

these problems, any system meant to perform this type of calculations must be 
able to manipulate algebraic and differential expressions. For this reason, such sys- 
tems are generally implemented as extensions, or packages, for a general purpose 
computer algebra system: the latter takes care of handling expressions, and the ex- 
tension introduces the differential-geometry specific features — differential forms, 
^ tensors, connections and so on. Examples are given by the packages difforms and 

GRTensor (see [9]) for Maple, or the Ricci package for Mathematica. 

A remarkable consequence of this approach is that one is essentially limited, when 

(T^ implementing one's own algorithms, to using the programming language built in the 

_il computer algebra system. Some drawbacks and limitations of these programming 

f^ languages are described in |T|, where an alternative is also introduced, namely the 

00 C++ library GiNaC. As suggested by the name (an acronym for GiNaC Is Not A 

^^ CAS), GiNaC differs from the above mentioned computer algebra systems in that 

^ it is based on a general purpose, well-established programming language such as 

C++, rather than introducing a new one. In particular, development, debugging 

^ and documentation of a program based on GiNaC can take advantage of the many 

^ tools commonly available to a C++ programmer. This makes GiNaC a natural 

choice when implementing new, complex algorithms, like the one introduced in [5], 

which provided the original motivation for the present work. 

In this paper we introduce Wedge, an extension to GiNaC that can be used to 
write a C++ program that performs computations in differential and Riemannian 
geometry. Wedge is able to perform algebraic or differential computations with 
differential forms and spinors, as well as curvature computations in an adapted 
frame, and contains some support for vector spaces, represented in terms of a basis. 
Bases of vector fields, or frames, play a central role in Wedge: in particular, the 
tangent space of a manifold is represented by a frame, and a Riemannian metric 
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2 DIEGO CONTI 

is represented by a (possibly different) orthonormal frame. Notice that the above- 
mentioned package GRTensor also supports working with adapted frames, but our 
approach differs in that the frame need not be defined in terms of coordinates; this 
can be useful when working on a Lie group, where the geometry is defined natu- 
rally by the structure constants, or on the generic manifold with a fixed geometric 
structure (see Section [2]). 

Another unique feature of Wedge among packages for differential and Riemann- 
ian geometry, beside the choice of C++, is the fact that it is completely based on 
free, open-source software. Like GiNaC, Wedge is licensed under the GNU general 
public license; its source code is available at http://libwedge.sourceforge.net,. 

This paper is written without assuming the reader is familiar with C++, and its 
purpose is twofold: to introduce the main functionality of Wedge, and to illustrate 
with examples certain features of C++ which can prove very helpful in the practice 
of writing a program to perform some specific computation. 

In the first section we introduce some basic functionality, concerning differential 
forms and connections; at the same time, we illustrate classes and inheritance. 

In the second section we introduce spinors and "generic" manifolds; at the same 
time, we illustrate object-oriented programming. 

In the third section we explain briefly how GiNaC handles expressions, then 
introduce bases and frames, namely the linear algebra features of Wedge. 

In the final section we give an application from Cartan-Kahler theory, with a 
short program that reproduces the computations of [3J that prove the local existence 
of metrics with holonomy G2 . 



1. Working with differential forms and connections 

A fundamental feature of C++ is the possibility of introducing user-defined types, 
i.e. classes (or structs). The definition of a class specifies not only the type of data 
contained in a variable which has that class as its type, but also some operations 
that can be performed on this data. This is generally better than having global 
functions which either take a long list of arguments or use global variables, and this 
can be seen, for instance, in situations where more sets of data appear in the same 
program. In this section we shall illustrate this point with an example, considering 
a problem concerning multiple connections on a fixed manifold. In the course of 
the section, we shall introduce some of the essential functionality of Wedge. 

The basic geometric entity in Wedge is the class Mcoiif old. A variable of type 
Manifold represents a manifold in the mathematical sense, which is assumed to be 
parallelizable, and represented by a global basis of one- forms e^, . . . , e", with dual 
basis of vector fields ei, . . . , e„. This assumption is tailored on Lie groups, but one 
can also think of a Mgoiif old object as representing a coordinate patch; also, every 
manifold can be written as a quotient MjG of a parallelizable manifold, so one 
can always reduce to the parallelizable case. For instance, one can compute the 
curvature of a Riemannian metric on MjG applying the O'Neill formula (see [TU]'). 

As an example, we shall consider the nilpotent Lie group AT, characterized by 
the existence of a global basis of one- forms e^, . . . , e^ such that 



(1) de^ = 0, de^ = 0, de' = e^^^ de" = 
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Figure 1. Inheritance graph for classes inheriting from Manifold 



In the above equation, e^^ stands for the wedge product e^ A e'^, and so on; this 
shorthand notation is also used by Wedge, and so it will appear throughout the pa- 
per in both formulae and program output. Equations (fTj) translate to the following 
code: 

struct X : public ConcreteManif old , public Has_dTable { 
X() : ConcreteManifold(4) { 

Declare_d(e(l) ,0) ; 

Declare_d(e(2) ,0) ; 

Declare_d(e(3) ,e(l)*e(2)) ; 

Declare_d(e(4),e(l)*e(3)); 
} 

>; 

The first line means that X inherits from the two classes ConcreteManif old 
and Has_dTable. Roughly speaking, this means that it inherits the functionality 
implemented by these two classes; inheritance defines a partial ordering relation, 
and is best represented by a graph (see Figure [T| . Specifically, inheriting from 
ConcreteManif old ensures that the forms e^, . . . , e^ are defined (where the dimen- 
sion 4 appears as a parameter on the second line). The fact that X inherits from 
Has_dTable means that the operator d on the manifold X is known in terms of its 
action on the basis e^, . . . , e*; the relations ([I]) are given with calls to the member 
function Declare_d. These calls appear in the body of the constructor of X, which 
is invoked automatically when a variable of type X is constructed. We can now 
instantiate a variable of type X and perform computations with it, e.g. 

X M; 

cout«M . d (M . e (4) ) «endl ; 

has the effect of printing e^^. Notice that both the forms e* and the operator d are 
implemented as members of X; this means that every variable of type X has its own 
set of forms e^, . . . , e" and d operator, and so the code must tell the compiler which 
X should be used. This is achieved here by the "M." appearing in the function calls; 
however, inheritance provides an equivalent, more appealing alternative, as in the 
following: 
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Struct SomeCalculations : X { 

SomeCalculationsO ■[cout<<d(e(4) )<<endl; } 
} M; 

The class Manifold also implements the Lie bracket and Lie derivative as mem- 
ber functions, whilst the exterior and interior product of forms are implemented 
by global functions. Covariant differentiation and curvature computations are ac- 
counted for by the class Connection, which we now introduce with an example. 

Every almost-complex manifold (Af , J) admits an almost complex connection V 
such that its torsion form Q and the Nihenjuis tensor are related by 

me{X,Y) ^ N{X,Y) ; 

by [S] V can be obtained from an arbitrary torsion-free connection V by 

(2) 

VxY = VxY ~ Q{X, Y), 4Q{X, Y) - {VjyJ)X + J{{VyJ)X) + 2J{{VxJ)Y). 

Suppose one wants to compute such a connection in the case of our nilpotent Lie 
group X, with the almost-complex structure determined by 

-'^(ei) = 62, J(e3) = 64. 

We can do so with the following code: 

struct AlmostComplex : public X ■[ 
TorsionFreeConnection<true> h; 

ex J(ex Y) {return Hook(Y,e(l) *e(2)+e(3)*e(4) ) ; } 
ex A(ex X,ex Y) { 

return h . Nabla<VectorField> (X , J (Y) ) - J (h . Nabla<VectorField> (X , Y) ) : 
} 
ex Q(ex X, ex Y) { 

return (A(J(Y) ,X)+J(A(Y,X) )+2*J(A(X,Y) ) )/4; 
} 

AlmostComplex : omega(this ,e() ) { 
Connection k(this,e()); 
for (int i=l;i<=4;++i) 
for (int j=l; j<=4;++j) 

k.DeclareNabla<VectorField>(e(i) ,e(j) , 

h.Nabla<VectorField>(e(i) ,e(j) )-q(e(i) ,e(j))) ; 
cout«k<<k.Torsion() ; 
} 
} M; 

This code defines and instantiates a struct named AlmostComplex inheriting from 
X; AlmostComplex has a data member h of type TorsionFreeConnection<true>, 
which represents a generic torsion- free connection V on X. The constructor uses 
h to compute another connection k, corresponding to V in Q, then prints the 
connection forms of V and its torsion form O. The latter is a vector 

(0, 0, -i632-ie«, \e^^~\e^') 

consistently with the fact that the almost complex structure J is not integrable. 

There are two connection objects appearing here, both represented internally by 
their connection forms in terms of the frame 6i, . . . , 6„. On construction, they are 
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initialized in terms of symbols Tijk, representing the functions 

where (, •) represents the pairing on TM (g) T*M . Since the type of h is 
TorsionFreeConnection<true>, the constructor of h ensures that the torsion is 
zero by solving the linear equations in the Tijk given by 

(3) ^e*AVe.e^-de^ =0. 

i 

Of course this does not determine V uniquely, and some of the Tijk remain as 
parameters. The conditions (pi) on the connection k are imposed directly by calling 
the member function DeclareNabla. The tensor Q is computed by the function Q, 
which in turn calls Nabla; since the latter is a member function, the connection h is 
used to compute the covariant derivative. Thus, the connection forms of k depends 
on the Fijfc, although the torsion Q is independent of the Tijk, as expected from 
the theory. 

Another important point is that DeclareNabla does not assign values directly to 
the connection forms. Rather, it uses GiNaC's function Isolve to solve a system 
of equations, linear in the Tijk, and substitutes the solution into the connection 
forms. This procedure is more flexible than giving the list of the Fyfe for two 
reasons: the arguments passed to DeclareNabla need not be elements of the frame 
(they can even be forms of degree higher than one, or, for Riemannian connections 
only, spinors), and part of the connection forms may be left unspecified. However, 
said procedure only works because the Tij^ are stored in the Connection object, 
which is how DeclareNabla knows it should leave the Fyfc as parameters, and solve 
with respect to the Tijk- 

Summing up, we have illustrated how the fact that a connection is represented 
as an instance of a class makes working with two different connections on a fixed 
manifold as natural as defining two variables of the same type. 

2. Torsion free connections and generic manifolds 

Class X from Section [l] represents a fixed manifold, where the action of the d 
operator can be recovered from its action on a basis of one-forms; we have seen that 
Wedge enables one to perform torsion computations, or impose torsion conditions. 
One can also go the other way, and compute the action of d in terms of the covariant 
derivative, with respect to a torsion-free connection. To illustrate this, suppose one 
has a four-dimensional Riemannian manifold X, with a global orthonormal frame 
e^, . . . , e^. The parallelism induces a spin structure; in general, on a parallelizable 
n-dimensional manifold a spinor can be viewed as a map X — > S, where E is the 
n-dimensional spinor representation. Let m — [n/2] be the integer part of n/2; E 
is a complex vector space of dimension 2™, represented in Wedge in terms of the 
basis uo, ■ • • , "2™-!, where u^ corresponds to 

ra — l 

u((-ir-\...,(-ir°), fc= ^(1-6^)2'- 

in the notation of [2] (see also [6] for more explicit formulae). In our case, we can 
declare that the constant spinor ij) : X ^ Yi, if) = uq \s parallel, i.e. 

Ve,^-0 z = l,...,4 
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with the following code: 

struct X : public Manif oldWith<RiemannianStructure> {. 
X() : ManifoldWith<RiemannianStructure>(4) { 
for (int i=l;i<=4;++i) 

DeclareNabla<Spinor>(e(l) ,u(0) ,0) ; 
> 

>; 

Hence, X is a generic parallelizable Rieniannian 4-nianifold with holononiy 
contained in SU(2); this condition does not determine the connection form 
uniquely, but it does impose certain conditions on them. Internally, the tem- 
plate class Manif oldWith<RiemajinianStructure> contains a member of type 
LeviCivitaConnection<f alse>, which behaves in a similar way to its counter- 
part TorsionFreeConnection<true> from Section IT] with two differences: first, it 
represents the Levi-Civita connection, which is both torsion-free and Riemannian, 
and secondly it does not use the operator d of the manifold to impose the condition 
(P = 0. Instead, Manif oldWith uses the Levi-Civita connection to compute the 
action of d via (l3| . 

Convention. All further code fragments appearing in this section will be as- 
sumed to appear in the body of the constructor of a class such as SomeCalculations 
of Section [TJ 

In fact, whilst the new X does not derive from has_dTable, the definition of 
SomeCalculations of Section[T|is still legitimate. Of course, the output will depend 
on the Tijk which in part are left unspecified. In particular, it is not possible to 
guarantee that d^ = 0, as this leads to an underdetermined system of differential 
equations, quadratic in the r^j. and linear in their derivatives. However, there 
are many interesting relations not involving the derivatives of the r^fc that can be 
proved using Manif oldWith. As a first example, we can verify that the existence of 
a parallel spinor is equivalent to the existence of a local frame e^, . . . , e* such that 
the forms 

(4) e'^ + e'\ e^^ + e^\ e^^ + e^\ 

are closed (see |7j). In our case, since we have chosen mq £^s the parallel spinor, the 
standard frame e^, . . . , e^ satisfies this condition, and the code 

cout«d(e (1) *e (2) +e (3) *e (4) ) «endl ; 
cout«d(e (1) *e (3) +e (4) *e (2) ) «endl ; 
cout«d(e (1) *e (4) +e (2) *e (3) ) «endl ; 

prints three times zero. One can obtain the opposite implication by invoking 
Declare_d in the constructor of X to impose that the forms Q are closed. 

Compared to the examples of Section [l] this code shows an important feature of 
C-H- which is common to many programming languages, but not so common among 
computer algebra systems: one can define different functions with the same name, 
and the compiler is responsible for selecting the correct one based on the context. 
In this case the functions d and Declare_d appear as members of different classes, 
i.e. Has_dTable and Manif oldWith, and so the meaning of a call to d in the body 
of a member of X depends on whether X inherits from Has_dTable, as was the case 
in Section [T] or Manif oldWith, as in the example above. In fact, C++ allows even 
greater flexibility, as we now illustrate with a second example. 
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In order to give a slightly more complicated application, we observe that our 
manifold X is in particular an almost-Kahler manifold with respect to any one of 
the closed two-forms Q; we shall fix the first one. Thus, it makes sense to consider 
the bilagrangian splitting TX = F (B G, where F = (ei, 63) and G = (e2, 64); this 
splitting determines a canonical connection w. By [11], the torsion of u; is zero if 
and only if the distributions F and G are integrable. We can actually prove this 
equivalence with Wedge; we shall illustrate the "only if" implication here. 

Connection omega (this ,e() , "Gcmima' ") ; 
for (int k=l;k<=4;++k) { 

omega. DeclareNabla<DifferentialForm>(e(k) ,e(l)*e(2)+e(3)*e(4) ,0) ; 

for (int i=l;i<=4;++i) 

for (int j=i7.2+l; j<=4; j+=2) 

omega.DeclareZero(Hook(e(j) , 

omega. Nabla<DifferentialForm>(e(k) ,e(i)))) ; 
} 

for (int k=l;k<=4;++k) 
for (int i=k7.2+l;i<=4;i+=2) 
for (int j=k7.2+l; j<=4; j+=2) 

omega. DeclareZero(Hook(e(j) , 

omega. Nabla<VectorField>(e(k) ,e(i) )-LieBracket (e(k) ,e(i) ) ) ) ; 
exvector T=omega.Torsion() ; 
DeclareZero(T.begin() ,T.end()) ; 
cout«LieBracket(e(l) ,e(3))«" , "«LieBracket (e(2) ,e(4))«endl; 

This code defines a generic connection lu, whose connection parameters are de- 
noted by r^^^ to distinguish them from the parameters Tijk of the Levi-Civita 
connection, and imposes the two sets of conditions that determine the symplectic 
connection. First, the holonomy is reduced to GL(2,M) by requiring that the sym- 
plectic form be parallel and the two Lagrangian distributions F and G preserved 
by the covariant derivative. Then one imposes 

Vx,>G = [^f,n?]G, '^Xa[YF]^[XG,YF]F, X,YeT{TX), 

where the subscripts denote projection (reflected in the code by the use of the in- 
terior product function Hook). Notice that LieBracket is a member function of 
Manifold, which Manif oldWith reimplements in terms of the Levi-Civita connec- 
tion; thus, its result depends on the Tijk- By the general theory, the conditions 
determine lu completely, i.e. lo does not depend on the F^ j, but only on the Fy^. 
Since this dependence is linear, the condition that the torsion is zero gives equa- 
tions in the Tijk that can be solved by the function DeclareZero, a member of 
Manif oldWith. Having imposed this torsion conditions, the program concludes 
that [61,63] and [62,64] equal 

eiFi42 + 63F331, F43ie4 + F23162 

respectively, proving that the distributions F and G are involutive. 

The code at work here uses a typical object-oriented programming technique. 
Specifically, Connection needs to compute the action of d in order to compute the 
torsion; to this end, every Connection object contains a pointer to the Manifold 
object it refers to. In this case, the manifold to which a; refers is represented 
by a Mcoiif oldWith object, which implements d using the Levi-Civita connection; 
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since d is a virtual function, the calls to d performed by Torsion execute the 
implementation of Manif oldWith. In the examples of Section [T] those same calls 
execute the implementation of Has_dTable. This ensures that the torsion of lu 
is computed correctly in terms of the Tijk here, whereas the code of Section [l] 
computes the torsion using the action of d on the forms e^, . . . ,e^, although the 
function Torsion is exactly the same. 

Remark. There is nothing essential about the assumption that X has holonomy 
SU(2). In fact, one can easily modify the code to obtain the same result for a 
generic symplectic 4-manifold X. 

3. Forms and frames 

In this section we introduce the main linear algebra functionality of Wedge; 
before doing that, we need to explain how GiNaC handles expressions. 

We have seen in Section [T] that manifolds are represented in Wedge by a global 
basis of one-forms e^, . . . , e", which can be accessed via a member function of the 
class Mcoiifold. The C-H- type of a form, say X.e(l), is the class ex defined in the 
library GiNaC, which handles expressions such as 2*e(l)+e(2), the result of which 
is still an ex. The class ex acts as a proxy to the GiNaC class basic. In practice, 
this means that an ex contains the address of an object whose type is a class 
that inherits from basic (e.g. add, representing a sum, or numeric representing a 
number), and it is this object that performs the actual computations. The use of 
virtual functions ensures that the correct code is used, depending on the C-H- type 
of the object the ex points to (for details, see pP). 

Thus, while X.e(l) has type ex, under the hood it "refers" to an object of type 
Dif f erentialDneForm, and it is this type that implements skew-commutativity, 
making use of a canonical ordering that exists among basic objects. More signifi- 
cantly, objects of type Dif f erentialDneForm and linear combinations of them can 
be used interchangeably; nonetheless, they can be distinguished by their internal 
representation, since the former have type Dif f erentialDneForm and the latter 
have type add (e.g. e(l)+e(2)) or mul (e.g. 2*e(l)). We refer to objects in the 
former category as simple elements. Mathematically, we are working with sparse 
representations of elements in the vector space generated by all the variables of 
type Dif f erentialDneForm appearing in a program. All this works equally well 
with other types than Dif f erentialDneForm. 

In general, given a small set of vectors, it is possible to switch from sparse to dense 
representation by extracting a basis, and writing the other terms in components. 
This is accomplished by the class template Basis. As a simple example, suppose X 
defines the Iwasawa manifold, characterized by the existence of an invariant basis 
e^, . . . ,e^ satisfying 

de^ ^ = de^ = de^ ^ de\ de^ = e^^ + e^^, de'' = e^^ + e^\ 

We can compute a basis of the space of invariant exact three-forms with the fol- 
lowing code (with the usual convention of p. O: 

Basis<Dif f erentialForm> b; 
for (int i=l;i<=6;++i) 
for (int j=i+l; j<=6;++j) 

b.push_back(d(e(i)*e(j) ) ) ; 
cout<<b; 
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resulting in the output 

g412^ _g312^ _g342^ g341^ g642 _ ^631 _ ^352 ^ ^415^ 

One can then recover the components of, say, de*^ by 

cout<<b. Components (d(e (4) *e (5) ) ) ; 

obtaining the vector ( o o o -i o ). 

Another use of Basis is related to the identification of each vector space with its 
dual. Indeed, the existence of simple elements enables us to establish a canonical 
pairing, which is bilinear and satisfies 




a,P) — i , a,P simple. 



In particular, we can identify vector fields with one-forms, by interpreting the pair- 
ing (,) as the interior product. For instance, the standard frame e^, . . . ,e" asso- 
ciated to a ConcreteManif old always consists of simple elements, and so in this 
case the form e* and the vector field e^ are represented by the same element e(i). 
However, this fact only holds for (,)-orthonormal bases. In general, the dual basis 
has to be computed, in order to deal with equations such as ([s]), where both a basis 
of one-forms and the dual basis of vector fields appear; the member dual of Basis 
takes care of this. 

Basis is based on the standard container vector, part of the Standard Template 
Library. So, a Basis object can be thought of as a sequence xi, . . . , a;„ of ex, where 
the Xk are linear combinations of simple elements of a given type T, which is a 
template parameter of Basis. However, Basis differs from vector in the following 
respects: 

• When elements are added to a Basis object, resulting in a sequence of 
possibly dependent generators xi, . . . , a;„, an elimination scheme is used to 
remove redunant (dependent) elements. This is done in such a way that 
the flag Vi C ... C Vn, 

Vk = Span{a;i,...,Xfc} 

is preserved. 

• The first time the member functions Components or dual are called. Basis 
sets up some internal data according to the following procedure. First, the 
set 

5* = S{xi, . . . , Xm) = {a simple | (a, Xk) 7^ for some 1 <k < m\ 

is computed; then, the sequence {xi, . . . , a;^} that represents the basis 
internally is enlarged with elements Xk £ S , m < k < n so that 

{xi, . . . , Xn} is a basis of Span S. 

Then, for each a G S* the components haj satisfying 

n 
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are computed; as a matrix, (baj) is the inverse of the matrix with entries 
Gja — {xj, a). The dual basis x"'^, . . . , x" is then given by 



X 

aes 



The actual code also take advantage of the fact that with a suitable ordering 
of the columns, one obtains the block form 

The dual basis and the baj are cached internally, until the basis is modified. 

• The member function Components uses the pairing (,) to determine the 
components with respect to S, then multiplies by the matrix (baj) to obtain 
the components with respect to the basis xi, . . . Xm- 

From the point of view of performance, the operations of Basis have the following 
complexity, measured in term of the number of multiplications: 

• Given a list Xi, . . . , a;„ of elements which are known to be linearly indepen- 
dent, one can construct a Basis object with no overhead over vector and 
no algebraic operation involved. 

• Inserting elements yi,...,yk in a Basis object consisting of elements 
xi, . . . ,Xn has a complexity of 0{m{n + k)'^), where 

m= \S{xi,...,Xn,yi,---,yk)\- 

• Setting up a Basis object consisting of elements xi, . . . , a;„ in order that 
the dual basis or components of a vector may be computed has a complexity 
of 0{{n + 171)71^), where 

m= |S'(a;i,...,a;„)|. 

• Once the Basis object is so set up, computing the components of a vector 
has complexity 0{nm). 

Since we work in the symbolic setting, the number of multiplications only gives 
a rough estimate of execution time. In any case, the above figures show the im- 
portance of being able to define a class that only performs certain operations when 
it is needed. Again, this is something which can be achieved naturally in C-H-, by 
storing a flag in each object Basis. Since a Basis object is only modifled by in- 
voking certain member functions, the flag is cleared automatically when the object 
is modified, effectively invalidating the cached values of baj. 

Remark. The frames associated to a Manifold, RiemannianStructure or 
Connection object are represented by Basis objects; thanks to the mechanisms 
described in this section, these frames need not consist of simple elements. In fact, 
this was the main motivation for introducing Basis. 

4. An APPLICATION: Cartan-Kahler theory 

In this section we consider, as an application, the problem of determining whether 
a linear exterior differential system (EDS) is involutive. This problem eventually 
boils down to computing the rank of certain matrices, but actually writing down 
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these matrices requires a system that supports differential forms and interior prod- 
ucts, as well as linear algebra; this makes Wedge particularly appropriate to the 
task. 

Suppose we have a real analytic EDS with independence condition (I, 0i, . . . , 9k) 
on a manifold X. As a special case, we shall take X to be the bundle of frames over 
a 7-manifold M, the di the components of the tautological form and T generated 
by the exterior derivative of the G2-forms (f) and *0. We wish to determine whether 
the system is involutive, i.e. every point of X is contained in some submanifold 
S C X such that all the forms of T vanish on S and the forms 0i\s, ■ ■ ■ ,0k\s are a 
basis of the cotangent bundle T*S; we then say that S is an integral manifold for 
2. In the special case, the EDS turns out to be involutive, proving the existence of 
local metrics with holonomy G2 (see [3 ). 

The general procedure is the following. Complete 61,..., 9^ to a basis 
01, . . . , 0fe, wi, . . . , Wn of 1-forms on X. Suppose that E C T^X is a fc-dimensional 
subspace such that the forms 9i\^ . . . ,9k are independent on E; then one can write 

We say that E is an integral element if the forms of T restrict to on E. The 
{x,Pij) are coordinates on the Grassmannian of fc-planes over X, and by definition 
the tangent space of an integral manifold at each point will be an integral element 
E of the form above. The integral elements of dimension n form a subset Vn of the 
Grassmannian. 

Suppose that there exists an integral element E C T^X, and Vn is smooth about 
E. By the Cartan-Kahler theory, the (local) existence of an integral manifold S 
with T^S = E is guaranteed if one can find a flag 

EqC ...C Ek = E, dimEk = k, 

such that 

dim K = Co + . . . + Cfe , 
where Cj is the dimension of the space of polar equations 

H{E,) = {a(ei, ...,€,,-) \ a e I, ei e E,} C T^X. 

Since we are assuming that Ej is contained in an integral element E, H{Ej) has 
the same dimension as its image H*{Ej) under the quotient map 

We shall say that a system is linear if 

X C A* Span{e^, . . . , e '} A Spanjcji, . . . ,Ci;„}, 

where the span is over C°°{X). This condition implies that the equations defining 
Vn are affine in the pij. In addition, one has H*{Ej) = H*{E'A, where _E' is the 
horizontal projection of Ej determined by the splitting 

{9i,...,9k,)®{uJl,...,UJn). 

One can compute the dimension of Vn and the Ci with the surprisingly short 
program of Fig. [2] This program represents the frame bundle P as a parallelizable 
manifold with a frame {9i,ujjk), where the 9i are the components of the tautological 
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form, and the LUjk are the components of a torsion-free connection form, related by 
the structure equation 

de' = eJ A ujij . 
The member function GetEquationsForVn simply replaces every occurrence of 
Ujk with pjk in terms of the coordinates Pij of the Grassmannian, and im- 
poses that the coefficients of the resulting forms are zero. The member function 
GetReducedPolarEquations, on the other hand, calls itself recursively to obtain all 
the forms e^^j . . .jei,ja, where a ranges in a basis of X and 1 < Ji < . . . < i; < j. 
When the resulting degree is one, the function stops the recursion and stores the 
differential one-form, applying a substitution corresponding to the projection ([5]). 

The code of Fig. [2] is for a general linear EDS on the bundle of frames. In order 
to obtain a result for the special case of G2, we can use the following: 
struct G2 : EDS { 

G2() : EDS (7) { 
1st I; 
I=d(ParseDifferentialForm(e(), "567-512-534-613-642-714-723")), 

dCParseDifferentialFormCeO, "1234-6712-6734-7513-7542-5614-5623")); 
Basis<symbol> A; 
GetEquationsForVn (A, I) ; 
for (int j=0;j<7;++j) 
{ 

Basis<Dif f erentialForm> V; 

for (1st : : const_iterator i=I .beginO ; i ! =1 . endO ;++i) 

GetReducedPolarEquations (V, *i , j) ; 
cout<<V. size()«endl; 
> 

cout<<A . size <<endl ; 
} 

}; 

The output shows that the Ci are 

0,0,0,1,5,15,28 

whereas the codimension of Vn is 49; thus, Cartan's test applies and the system is 
involutive. 

Notice that the equations on the Grassmannian are stored in a Basis<symbol> 
object, and so the codimension of V7 is just the size of the basis. This makes 
sense because we know that the equations defining V7 are linear in the pij] how- 
ever, the container A is passed as a template parameter to GetEquationsForVn. 
This means that if we were dealing with an EDS where the equation are affine in 
the Pij rather than linear, e.g. the EDS associated to nearly-Kahler structures in 
six dimensions, it would suffice to modify the code by declaring A to be of type 
Af f ineBasis<symbol>. The general philosophy, used extensively in Wedge, is that 
of writing "generic" template functions, without making specific assumptions on 
the nature of the arguments, so that the type of the arguments determines the 
precise behaviour of the code. 

Remark. Concerning the non-linear case. Wedge also provides a class template 
PolyBasis, analogous to Basis, where polynomial equations can be stored, relying 
on CoOaA for the necessary Grobner basis computations (see |1]). Notice however 
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Struct EDS : public ConcreteManif old, public Has_dTable { 
1st moduloIC; 
int dim; 

EDSCint dimension) : ConcreteManif old(dimension*(dimension+l) ) { 
dim=dimension; 
for (int i=l;i<=dim;++i) { 
ex x; 
for (int j=l; j<=dim;++j) 

x+=e ( j ) *e (i*dim+j ) ; 
Declare_d(e(i) ,x) ; 
moduloIC. append(e(i)==0) ; 
> 
} 

template<typencmie Container> 

void GetEquationsForVn(Container& container, 1st 1) { 
1st substitutions; 

for (unsigned i=dim+l; i<=dim*(dim+l) ;++i) { 
ex ei ; 
for (unsigned j=l ; j<=dim;++j) 

ei+=symbol("p"+ToString(i,j))*e(j) ; 
substitutions . append(e(i)==ei) ; 
} 

for (1st : : const_iterator i=l.begin(); i !=1 . endO ;++i) 
GetCoefficients (container , i->subs (substitutions) ) ; 
> 

template<typen£mie Container> 

void GetReducedPolarEquations(Container& container, ex form, int j) { 
if (degree<Dif f erentialForm>(f orm)==l) 

container .push_back(form. subs (moduloIC) ) ; 
else if (j>0) { 

GetReducedPolarEquations (container , form, j-1) ; 
GetReducedPolarEquations (container , Hook(e(j) ,form) , j-1) ; 
} 



>; 



} 



Figure 2. Using Wedge to determine if a linear EDS on the frame 
bundle is involutive 

that the simple program of Fig. [2] assumes linearity when computing the polar 
equations. 
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